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ANALYTICAL STUDY OF MECHANISMS FOR 
NITRIC OXIDE FORMATION DURING COMBUSTION OF 
METHANE IN A JET -STIRRED COMBUSTOR 

Casimir J. Jachimowski 
Langley Research Center 

SUMMARY 

The role of chemical kinetics in the formation of nitric oxide during the combustion 
of methane has been examined analytically by means of a detailed chemical mechanism for 
the oxidation of methane, for the reaction between hydrocarbon fragments, and for the for- 
mation of nitric oxide. By comparing predicted nitric oxide levels with values reported in 
the literature from jet-stirred combustor experiments, it was determined that the nitric 
oxide levels observed in fuel-rich flames cannot be described by a mechanism in which 
the rate of nitric oxide formation is controlled solely by the kinetics of oxygen atom for- 
mation. A proposed mechanism for the formation of nitric oxide in methane-rich flames 
reproduces the observed levels. The oxidation of HCN appears to be an important factor 
in nitric oxide formation. 


INTRODUCTION 

Much attention has been paid recently to the development of analytical models that 
can predict quantitatively the amounts of nitric oxide emitted from combustion systems 
such as aircraft gas turbines. In most of these models, it is assumed that nitric oxide is 
formed primarily through the reactions 

(1) O + N2 - NO + N 

(2) N + 02 - NO + O 

(3) N+OH-NO + H 

with the rate of nitric oxide formation being limited by the rate of reaction (1). The rate 
of reaction (1) is determined by the flame temperature and the level of oxygen atoms which, 
in turn, are determined by the kinetics of the fuel-oxygen reactions. This reaction scheme 
for nitric oxide formation is often called the Zeldovich mechanism. 



In general, the current kinetic models are not able to adequately predict the observed 
nitric oxide levels, particularly for stoichiometric and fuel-rich hydrocarbon -air mixtures. 
For example, during an investigation of the formation of nitric oxide in hydrogen, carbon 
monoxide, and propane flames within a jet-stirred combustor, Engleman and co-workers 
(ref. 1) compared experimental nitric oxide concentrations with those predicted from a 
kinetic model of the appropriate system. Formation of nitric oxide in the nonhydro- 
carbon (H 2 and CO) systems was found to be adequately modeled at all equivalence ratios 
by reactions (1) to (3) and the appropriate fuel-oxygen kinetic reactions. However, the 
predicted levels in the propane flames greatly underestimated the observed concentra- 
tions (factors of 4 in fuel-lean and of 10 in fuel-rich flames). Similarly, Edleman and 
co-workers (ref. 2) compared experimental nitric oxide concentrations from methane-air 
flames in a jet-stirred combustor with those predicted from a kinetic model and found 
that the predicted levels underestimated the observed concentrations for stoichiometric 
and fuel-rich mixtures. It has been suggested (refs. 3, 4, and 5) that excess amounts 
C'superequilibrium") of oxygen atoms produced during the combustion of hydrocarbon fuels 
are responsible for the higher nitric oxide levels which are observed, and the current 
schemes do not adequately predict these oxygen atom levels. However, recent studies by 
Fenimore (ref. 6) and others (refs. 7 and 8) suggest that nitric oxide is produced by an 
additional reaction mechanism that only occurs in hydrocarbon -air flames. This nitric 
oxide is formed during the initial stages of the combustion reaction and is sometimes 
called "prompt" NO. It has been proposed that the reaction mechanism during early 
stages of combustion involves a hydrocarbon fragment which attacks molecular nitrogen, 
for example 


(4) CH + N 2 - HCN + N 

Nitric oxide could then be produced by reaction (2) which is very fast compared with 
reaction (1). 

The correct explanation for the discrepancy between predicted and observed nitric 
oxide levels in hydrocarbon flames is still a matter of controversy. However, it seems 
that a step toward resolving this controversy could be made if nitric oxide predictions are 
made by means of a hydrocarbon reaction model that reliably describes in a quantitative 
manner the combustion characteristics of the fuel, including the kinetics of oxygen atom 
formation. Comparison of these predictions with experimental data could then determine 
whether reactions (1) to (3) can adequately describe the nitric oxide formation process or 
whether an additional reaction mechanism must be used. Recently a study was made of 
the kinetics of oxygen atom formation during the oxidation of methane. (See ref. 9.) This 
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study resulted in a chemical kinetic scheme for methane oxidation which reasonably repro- 
duced the reported combustion characteristics of methane -oxygen mixtures behind shock 
waves, including the kinetics of oxygen atom formation. 

In this report the results of an analytical study of nitric oxide formation in methane - 
air flames within a jet-stirred combustor are presented. The purpose of the study was to 
compare predicted nitric oxide levels with experimental values by using the methane oxi- 
dation model reported in reference 9 coupled with reactions (1) to (3) and with a reaction 
mechanism which includes reaction (4) . 


SYMBOLS 

enthalpy of species i, kJ/mol (kcal/mol) 

AH-p heat of reaction at temperature T, kJ/mol (kcal/mol) 

k rate coefficient, cm^/mol-sec for bimolecular reactions, cm^/mol^-sec for 

termolecular reactions 

ih mass flow rate, g/sec 

N total number of species 

p pressure, atm (1 atm = 101.3 kPa) 

Q net heat loss from reactor, kilojoules (kilocalories) per gram of mixture 

T absolute temperature, K 

V reactor volume, cm^ 

Wj net rate of formation of species i, mol/cm^-sec 

cTj^ concentration, moles of species i per gram of mixture 


equivalence ratio. 


Fuel -air ratio 

Stoichiometric fuel -air ratio 


[N0],[HCN] concentration of nitric oxide and hydrogen cyanide, respectively, ppm 
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METHANE -AIR CHEMICAL KINETIC MODEL 


The chemical kinetic reaction model and the rate coefficients used in this study are 
listed in table I. The rate coefficients selected for these reactions were obtained from 
the literature whenever possible, (See refs, 9 to 25,) Rate coefficients that could not be 
found in the literature were estimated. Reactions (1) to (3) describe the formation of 
nitric oxide, and reaction (4) describes the formation of atomic nitrogen as suggested by 
Fenimore in reference 6, The rate coefficient for reaction (4) was estimated by using 
Hirschf elder's rule (ref, 26) to determine the activation energy, A preexponential factor 
of 10^^ cm^/mol-sec, common to most radical -molecule reactions, was selected for re- 
action (4), Reactions (5) to (23) describe the oxidation of methane. These reactions and 
rate coefficients provide a quantitative description of the kinetics of oxygen atom forma- 
tion and the ignition characteristics of methane -oxygen mixtures. Further details on this 
reaction scheme are given in reference 9, 

Reactions (25) to (30) describe the formation of the hydrocarbon fragments CH 2 and 
CH by the systematic removal of hydrogen atoms by the radicals H, O, and OH, The rate 
coefficients for reactions (25) to (27) were estimated by using activation energies reported 
for similar reactions (reactions (6) to (8)) and preexponential factors that are smaller by 
a factor of about 10, These preexponential factors were chosen so that the rate of these 
reactions would not be larger than the rate of reaction (10) which together with reaction (9) 
is considered to be the primary process responsible for the destruction of methyl radi- 
cals. (See ref. 9.) Reactions (31) to (33) describe the destruction of the hydrocarbon frag- 
ments CH 2 and CH by reaction with molecular oxygen. Reaction (33) is very exothermic 
(aH 2000 K ~ "356 kJ/mol (-85 kcal/mol)); therefore the activation energy was assumed to 
be zero. Since reaction (33) is similar to reaction (4), the same preexponential factor 
was chosen. Other reactions involving hydrocarbon fragments and molecular nitrogen 
are possible (ref. 27); however, they are highly endothermic and therefore would not be 
expected to be important compared with reaction (4). 

JET-STIRRED COMBUSTOR MODEL 

An ideal jet-stirred combustor consists of a combustion chamber with a well-defined 
volume into which the fuel -air mixture enters and instantaneously mixes with the reactor 
contents at a constant pressure. The mixture of reactants and products flows out contin- 
uously in a manner which results in a steady -state operation. In steady -state operation 
the species concentration within the reacting volume (combustion chamber) is determined 
from the balance between the net rate of production (or consumption) of each species by 
chemical reaction and the difference between the input and output flow rates of the species. 
This type of combustor is called a perfectly stirred reactor. A perfectly stirred reactor 



is treated mathematically by formulating the equations for mass, species, and energy con- 
servation as described by Jones and Prothero (ref. 28) in solving a system of simultaneous 
nonlinear algebraic equations. 

Briefly, the conservation equations for a perfectly stirred reactor are formulated 
in the following manner. Consider a stirred reactor of volume V containing combustion 
gases at a pressure p and a uniform temperature T. For steady-state operation, con- 
servation of mass is given by 

rii* = ih 

where m* is the mass flow rate into the reactor and rh is the mass flow rate out of 
the reactor. The superscript _* always refers to incoming reactants. Conservation of 
species is given by 

^(ai-of)=Wi (i= 1, 2, 3, . . .,N) 

where ct* and are the input and output concentrations of species i, N is the total 
number of species, and is the rate of formation of species i per unit volume be- 

cause of chemical reactions. Conservation of energy is given by 

N N 

I <H* - ^ .TjH, = Q 

i=l i=l 

where is the enthalpy of species i and Q is the heat loss from the reactor per 
unit mass because of radiation, thermal diffusion, and convection. For adiabatic opera- 
tion, of course, Q is zero. 

The stirred-reactor calculations presented in this report were made with a com- 
puter program described in reference 29. These calculations simulated the experiments 
made by Bartok and co-workers (ref. 30), who carried out experiments in a Longwell- 
Weiss jet-stirred combustor for fuel-lean, stoichiometric, and fuel-rich methane-air 
mixtures. Inlet temperature for these experiments was 463 K and reaction pressure was 
1 atm. The nitric oxide data selected for comparison were from experiments with resi- 
dence times of 2.3 msec. These data were selected because the residence times were 
short enough to increase chances that the reactor was nearly perfectly stirred, and also 
because a wide range of inlet equivalence ratio was used. Nitric oxide data were also 
reported for 1.8 -msec residence times; however, the equivalence -ratio range was not as 
wide as the 2.3 -msec data. 


5 


I 



DISCUSSION OF RESULTS 



The results presented in figure 1 are a comparison between the two processes 
describing nitric oxide formation. When the principal methane oxidation mechanism, 
reactions (5) to (24), is coupled with reactions (1) to (3), the calculated nitric oxide levels 
are in good agreement with the fuel-lean data only. For fuel-rich conditions, the pre- 
dicted levels are up to a factor of 10 lower than the observed levels. Parametric studies 
with the methane oxidation scheme revealed that reactions (5), (9), and (10) had a large 
influence on the predicted oxygen atom levels, particularly the peak oxygen atom levels. 
(See ref. 9.) An attempt was made to adjust some of the rate coefficients for these re- 
actions to provide better agreement between experimental and calculated results for the 
fuel -rich regime. However, a set of rate coefficients could not be found that would sub- 
stantially increase the oxygen atom levels. It seems that some other chemical path is 
required for fuel -rich conditions. 

When reaction (4) and reactions (25) to (33) are added to the principal methane oxi- 
dation scheme, the predicted nitric oxide levels are brought into fair agreement with the 
data, as shown in figure 1. However, the predicted nitric oxide levels are still smaller 



10 1 L 1 -i 1 

0.6 0.8 l.q 1.2 1.4 1.6 

Equivalence ratio, <p 


Figure 1.- Calculated and experimental nitric oxide levels for methane -air mixtures. 

"^initial ~ “^63 K; p = 1 atm. 
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by a factor of about 2 for 0 = 1,4. The addition of these reactions did not substantially 
change the nitric oxide levels for the fuel -lean and near stoichiometric regime because the 
rates of reactions (31) to (33) are sufficiently fast to prevent the buildup of the CH 2 and CH 
hydrocarbon fragments, and thus the rate of nitrogen atom production by reaction (4) is 
reduced. An attempt was made to adjust the rate coefficients for the fragment -forming 
reactions and the fragment -consuming reactions to provide better agreement between the 
experimental and calculated nitric oxide levels. However, a set of rate coefficients could 
not be found that would significantly increase nitric oxide levels for the very fuel-rich 
regime (<^ » 1.4) without overpredicting the nitric oxide levels in the near stoichiometric 
region. For example, by decreasing the rate coefficient for the fragment -consuming re- 
action by a factor of 10, there is better agreement between experiment and prediction at 
^ « 1.4; but at (f> = 1.1, the predicted value is about 50 percent larger than the measured 
level. This comparison is given in figure 2. It appears that the kinetic model presented 
in table I is not complete and other sources of nitric oxide must be found. 

An additional source of nitric oxide could be the hydrogen cyanide, HCN, that is 
produced in reaction (4) . The kinetic model in table I does not provide any path by which 


400r- 


^(^Experimental data, ref. 30 



10! ^ ^ I 

0.6 0.8 1.0 1.2 1.4 1.6 


Equivalence ratio, <J) 


Figure 2.- Experimental and calculated nitric oxide levels using ^prompt” 
NO reactions in table I and "prompt^' NO reactions with adjusted rate 
coefficients for reactions (31) to (33). 


HCN can be consumed and permits the HCN to accumulate. The accumulated HCN levels 
can be substantial and approach the nitric oxide levels as shown in figure 3. It is likely 
that significant amounts of HCN can be oxidized to nitric oxide, particularly in a stirred 
combustor where products are continually being mixed with fresh reactants. If it is 
assumed that all the HCN can be oxidized to nitric oxide, the predicted nitric oxide levels 
are in excellent agreement with the observed nitric oxide levels, as shown in figure 4. 

The conversion of hydrogen cyanide to nitric oxide could take place directly through a 
reaction such as 

(34) HCN + O 2 * HCO + NO 

which is exothermic (AH 200 O K “63 kJ/mol (-15 kcal/mol)) or through an indirect process 
leading to the formation of CN which then reacts directly with molecular oxygen by the very 
exothermic reaction (AH 20 OO K = kJ/mol (-114 kcal/mol)) 

(35) CN+O 2 -NO + CO 



0.6 0.8 1.0 1.2 1.4 1.6 

Equivalence ratio, 4> 


Figure 3.- Calculated NO and HCN levels for methane-air mixtures. 
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400r- 


/\ Experimental data, ref. 30 



Figure 4.- Calculated NO and NO + HCN levels for methane -air mixtures. 


It is not claimed that the reaction mechanism presented for nitric oxide formation 
is unique or corresponds to reality. However, some reaction process other than the 
direct formation of nitric oxide by reaction (1) appears necessary for fuel-rich conditions. 
The hypothesized kinetic mechanism contains reactions that probably do occur to some 
extent, and the ability of this set of reactions to reasonably reproduce the experimental 
observations requires that the assigned rate coefficients be reasonably correct. Nonethe- 
less, experimental verification of the individual rate coefficients is still required. Also 
more information is required on the kinetics of HCN and CN oxidation. Hydrogen cyanide 
has been observed in hydrocarbon-air flames as an intermediate species and can exceed 
the nitric oxide levels in very fuel-rich flames. (See ref. 1.) Comparison of predicted 
HCN levels with levels measured in various stirred-reactor experiments with fuel-rich 
methane -air mixtures could proyide a means to verify the reaction scheme proposed in 
this report. Other hydrocarbon-air mixtures should also be studied to determine whether 
other reactions can contribute to nitric oxide formation in fuel -rich mixtures. Different 
reactions may become important because the chemical structure of the fuel is different. 
For example, combustion of unsaturated hydrocarbons like ethylene, acetylene, and aro- 
matics could result in larger HCN levels from reactions such as 
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(36) C2H + N2 ~ HCN + CN 

(37) C 2 + N 2 - 2CN 

As other sources of hydrocarbons, such as shale oil and various coal syncrudes, begin to 
be used in the production of hydrocarbon fuels, including jet fuels, the chemical composi- 
tion of available fuels will vary even more. Any attempt to reliably predict nitric oxide 
levels as a function of chemical composition will require a greater understanding of the 
effect of fuel structure on nitric oxide production, particularly in fuel-rich mixtures. 

CONCLUDING REMARKS 

The purpose of this study was to examine the role of chemical kinetics in the for- 
mation of nitric oxide during the combustion of methane in a jet-stirred combustor. The 
results of the study revealed that nitric oxide formation in fuel-lean methane -air mixtures 
can be explained by a kinetic model in which the rate of nitric oxide formation is controlled 
by the kinetics of oxygen, atom formation. However, for fuel -rich combustion another path 
for nitric oxide formation appears to be required. Fuel-rich predictions can be brought 
into fair agreement with the observed nitric oxide levels by assuming a reaction sequence 
that produces nitrogen atoms through a reaction between the hydrocarbon fragment CH and 
molecular nitrogen. Better agreement was obtained if the intermediate HCN was allowed 
to react completely to nitric oxide. The proposed reaction scheme is not necessarily 
unique, and therefore, further experimental studies, such as measurement of HCN levels 
in a jet-stirred combustor, should be made to verify the kinetic scheme. Also, studies 
of the kinetics of HCN oxidation are needed to determine the mechanism of HCN consump- 
tion. Studies of NO and HCN formation in hydrocarbons other than methane are also 
recommended to help determine whether other reaction paths can also lead to nitric oxide 
formation. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, Va. 23665 
October 23, 1975 
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TABLE I.- METHANE-AIR CHEMICAL KINETIC MODEL 


Reaction 

Rate coefficient, k 

(a) 

Ref. 

(1) O + Ng - NO + N 

1.36 X 10^^ exp(-37946/T) 

10 

(2) N + Og - NO + O 

6.40 X lO^T exp(-3145/T) 

10 

(3) N + OH - NO + H 

4.00 X 10l3 

11 

(4) CH + Ng - HCN + N 

1.00 X 10^3 exp(-8052/T) 

Estimate 

(5) CH 4 + M - CHg + H + M 

4.00 X lO^"^ exp(-44500/T) 

9 

( 6 ) H + CH 4 - CHg + Hg 

1.26 X 10^^ exp(-5989/T) 

12 

(7) OH + CH 4 - CHg + HgO 

3.00 X 10^3 exp(-3020/T) 

13 

( 8 ) 0 + CH 4 - CHg + OH 

2.00 X 10^3 exp(-4630/T) 

14 

(9) CHg + Og - CHgO + OH 

1.70 X 10 l 3 exp(-745/T) 

9 

(10) CHg + 0 - CHgO + H 

1.3 X 10^^ exp(-1006/T) 

13 

(11) OH + CHgO - HCO + HgO 

2.30 X 10^3 

13 

(12) OH + HCO - CO + HgO 

1.00 X 10^3 

9 

(13) HCO + M- H + CO + M 

2.00 X 10 ^ 3 ^ 0.5 exp(-14400/T) 

15 

(14) OH + CO - COg + H 

4.00 X 10 l 3 exp(-4026/T) 

16 

(15) CO + Og - COg + 0 

1.60 X 10^3 exp(-20634/T) 

17 

(16) 0 + Hg - OH + H 

2.07 X 10l4 exp(-6920/T) 

18 

(17) OH + Hg - HgO + H 

5.2 X 10^3 exp(-3271/T) 

16 

( 18 ) H + Og - OH + 0 

1.22 X exp(-8369/T) 

18 

(19) OH + OH - HgO + 0 

5.5 X 10^3 exp(-3523/T) 

16 

(20) H + H + M-Hg + M 

1.00 X lolS 

19 

( 21 ) H + OH + M - HgO + M 

8.4 X io21t"30 

20 

(22) CO + 0 + M - COg + M 

6.0 X 10^3 

21 

(23) Og + M- O + O + M 

1.32 X IO^It'^-'^ exp(-59400/T) 

22 

(24) Ng+M-N+N + M 

4.98 X iol8x-0-75 exp(-113200/T) 

22 

(25) H + CHg - CHg + Hg 

1.00 X 10^3 exp(-6039/T) 

Estimate 

(26) 0 + CHg - CHg + OH 

1.00 X 10 I 2 exp(-5032/T) 

Estimate 

(27) OH + CHg - CHg + HgO 

1.00 X 10^3 exp(-3020/T) 

Estimate 

(28) H + CHg - CH + Hg 

2.7 X lO^^T®-®’ exp(-12934/T) 

23 

(29) 0 + CHg - CH + OH 

1.9 X exp(-12934/T) 

23 

(30) OH + CHg - CH + HgO 

Obtained from k 3 Q = k 2 g 

24 

(31) Og + CHg - HCO + OH 

1.0 X 10l4 exp(-1862/T) 

25 

(32) Og + CH - CO + OH 

Obtained from k 32 = 0 . 4 k 2 g 

24 

(33) Og + CH - HCO + 0 

1.00 X 10^^ 

Estimate 


^The units for k are cm^/mol-sec for bimolecular reactions and cm®/naol^-sec for ter- 
molecular reactions. Activation energies are expressed in kelvins. Subscripts to k correspond 
to reaction numbers. 
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